Rubber friction: role of the flash temperature 
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When a rubber block is sliding on a hard rough substrate, the substrate asperities will exert 
time-dependent deformations of the rubber surface resulting in viscoelastic energy dissipation in the 
rubber, which gives a contribution to the sliding friction. Most surfaces of solids have roughness 
on many different length scales, and when calculating the friction force it is necessary to include 
the viscoelastic deformations on all length scales. The energy dissipation will result in local heating 
of the rubber. Since the viscoelastic properties of rubber-like materials are extremely strongly 
temperature dependent, it is necessary to include the local temperature increase in the analysis. At 
very low sliding velocity the temperature increase is negligible because of heat diffusion, but already 
for velocities of order 10~ 2 m/s the local heating may be very important. Here I study the influence 
of the local heating on the rubber friction, and I show that in a typical case the temperature increase 
results in a decrease in rubber friction with increasing sliding velocity for v > 0.01 m/s. This may 
result in stick-slip instabilities, and is of crucial importance in many practical applications, e.g., for 
the tire-road friction, and in particular for ABS-breaking systems. 
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1 Introduction 

Rubber friction is of extreme practical importance, 
e.g., in the context of tires, wiper blades, conveyor belts 
and seaisp). Rubber friction on smooth substrates, 
e.g., on smooth glass surfaces, has two contributions, 
namely an adhesive (surface) and a hysteretic (bulk) 
contribution^, The adhesive contribution results 

from the attractive binding forces between the rubber 
surface and the substrate. Surface forces are often domi- 
nated by weak attractive van der Waals interactions. For 
very smooth substrates, because of the low elastic moduli 
of rubber-like materials, even when the applied squeezing 
force is very gentle this weak attraction may result in a 
nearly complete contact at the interface 0, 0, 0, leading 
to the large sliding friction force usually observed^. For 
rough surfaces, on the other hand, the adhesive contri- 
bution to rubber friction will be much smaller because of 
the small contact area. The actual contact area between 
a tire and the road surface, for example, is typically only 
~ 1% of the nominal footprint contact area|jj |s 111 • Un- 
der these conditions the bulk (hysteretic) friction mech- 
anism is believed to prevail H Q. For example, the 
exquisite sensitivity of tire-road friction to temperature 
just reflects the strong temperature dependence of the 
viscoelastic bulk properties of rubber. 

When a rubber block is slid on a hard rough substrate 
the surface asperities of the substrate will exert fluctuat- 
ing forces on the rubber surface which, because of the in- 
ternal friction of the rubber, will result in energy transfer 
from the translational motion of the block into the irreg- 
ular thermal motion. This will result in a contribution 
to the friction force acting on the rubber block. The en- 
ergy dissipation will result in local heating of the rubber. 
Since the viscoelastic properties of rubber-like materi- 
als are extremely strongly temperature dependent, it is 
necessary to include the local temperature increase when 
calculating the friction force. At very low sliding veloc- 
ity the temperature increase is negligible because of heat 



diffusion, but already for velocities of order 10 -2 m/s the 
local heating may become very important. Here I study 
the influence of the local heating on the rubber friction. I 
show that in a typical case the temperature increase will 
result in a friction which decreases with increasing sliding 
velocity for v > 0.01 m/s. This may result in stick-slip 
instabilities, and is of crucial importance in many prac- 
tical applications, e.g., for the tire-road friction, and in 
particular for ABS-breaking systems. 

2 Surface roughness and macro-asperity contact 

The influence of surface roughness on rubber friction is 
mainly determined by the surface roughness power spec- 
trum C(q) defined by 



C(q) 
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iq-x 
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Here /i(x) is the substrate height measured from the av- 
erage plane defined so that (h) = 0. The (. . .) stand for 
ensemble averaging, or averaging over the surface area. 
We have assumed that the statistical properties of the 
substrate are translational invariant and isotropic so that 
C(q) only depend on the magnitude q — |q| of the wave 
vector q. Note that from (1) follows 
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so that the root-mean-square (rms) roughness amplitude 
ho = (/i 2 ) 1 / 2 is determined by 



(h 2 ) = J d 2 qC(q) = 27T^ 



dq qC(q) 



(2) 



In reality, there will always be an upper, qi, and a lower, 
qL, limit to the q-integral in (2). Thus, the largest pos- 
sible wave vector will be of order 2n/a, where a is some 
lattice constant, and the smallest possible wave vector is 
of order 2ir/L where L is the linear size of the surface. 
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FIG. 1: Surface roughness power spectrum of a surface which 
is self affine fractal for gi > q > go- The long-distance roll- 
off wave vector go and the short distance cut-off wave vector 
gi depend on the system under consideration. The slope of 
the logC — logg relation for g > go determines the fractal 
exponent of the surface. The lateral size L of the surface (or 
of the studied surface region) determines the smallest possible 
wave vector g_L = 2n/L. 



Many surfaces tend to be nearly self-affine fractal. A 
self- affine fractal surface has the property that if part 
of the surface is magnified, with a magnification which 
in general is appropriately different in the perpendicular 
direction to the surface as compared to the lateral di- 
rections, then the surface "looks the same", i.e., the sta- 
tistical properties of the surface are invariant under the 
scale transformation. For a self-affine surface the power 
spectrum has the power-law behavior 



C(q) 



-2(H+1) 



where the Hurst exponent H is related to the fractal di- 
mension Df of the surface via H = 3 — Df. Of course, 
for real surfaces this relation only holds in some finite 
wave vector region qo < q < qi, and in a typical case 
C(q) has the form shown in Fig. ^ Note that in many 
cases there is a roll-off wavevector go below which C(q) 
is approximately constant. 

Asphalt and concrete road pavements have nearly per- 
fect self-affine fractal power spectra, with very well- 
defined roll-off wave vector q n — 2ir/\ of order 
1000 m _1 , corresponding to Ao ~ 1 cm, which reflect 
the largest stone particles used in the asphalt. This is 
illustrated in Fig. [21 for two different asphalt pavements. 
From the slope of the curves for q > go one can deduce 
the fractal dimension Df w 2.2, which is typical for as- 
phalt and concrete road surfaces. 

Assume that an elastic solid with a flat surface is 
squeezed against a hard, randomly rough substrate. Fig. 
[3]shows the contact between two solids at increasing mag- 
nification £. At low magnification (( 1) it looks as 
if complete contact occurs between the solids at many 
macro- asperity contact regions, but when the magnifi- 
cation is increased smaller length scale roughness is de- 
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FIG. 2: The surface roughness power spectra C(q) for two 
asphalt road surfaces. 



tected, and it is observed that only partial contact oc- 
curs at the asperities. In fact, if there would be no short 
distance cut-off the true contact area would vanish. In 
reality, however, a short distance cut-off will always exist 
since the shortest possible length is an atomic distance. 

The magnification £ refers to some (arbitrarily) chosen 
reference length scale. This could be, e.g., the lateral size 
L of the nominal contact area in which case £ = L/X, 
where A is the shortest roughness wavelength components 
which can be resolved at magnification £. In this paper 
we will consider surfaces with power spectra's of the form 
shown in Fig. ^ and we will use the roll-off wavelength 
Ao = 27r/<7o as the reference length so that £ = Ao/A = 

q/qo- 

I now explain the concept of macro- asperity contact 
area, which is important for my treatment of the influ- 
ence of the flash-temperature on rubber friction. Assume 
that the surface roughness has the qualitative form shown 
in Fig. ^with a roll-off wave vector qo corresponding to 
the magnification £ = 1. In this case the macro-asperity 
contact is the contact region between the solids when the 
system is studied at low magnification £ = £ m w 2 — 5 
(see Appendix A). At this magnification, at low squeezing 
pressures one observe a dilute distribution of randomly 
distributed macro-asperity contact regions with lateral 
size typically of order ~ Ao/Cm- As the magnification is 
increased each macro-asperity contact region breaks up 
into a relative dense distribution of much smaller micro- 
asperity contact regions. This is illustrated in Fig. 0] 
which shows the result of a Molecular Dynamics (MD) 
calculation 10] where the surface roughness power spec- 
trum was assumed to be of the form shown in Fig. ^ 
(with the fractal dimension 2.2). 

Consider the contact between two solids at low nom- 
inal contact pressure o~q = F^/Aq, where the contact 
area is proportional to the load. Consider first the sys- 
tem on the length scale Ao/Cm- On this length scale the 



solids will make (apparent) contact at a low concentra- 
tion of (widely separated) contact areas. Since the sep- 
aration between these macro-asperity contact regions is 
very large we can neglect the interaction between the 
macro contact regions: in this case the pressure in the 
macro-asperity contact regions will be of order ~ q^hoE, 
where ho is the rms roughness amplitude, and E the 
elastic modulus. Thus, the (average) pressure in the 
macroasperity contact regions is independent of the nom- 
inal contact pressure cjq = F^/Aq. Now, each macro- 
asperity is covered by smaller micro asperities, and the 
smaller asperities by even smaller asperities, and so on. 
It is easy to see that at short enough length scale the 
micro-asperity contact regions will be very closely sepa- 
rated and it is therefore impossible to neglect the (elastic 
or thermal) interaction between the asperities at short 
enough length scale. 

In Sec. 4 we will calculate the flash temperature in the 
asperity contact regions as a rubber block is sliding on 
a rough substrate. We will neglect the thermal coupling 
(or overlap) between the macro-asperity contact regions, 
which should be a good approximation as long as the 
contact area at the magnification C ~ Cm is much smaller 
than the nominal contact area (which implies that the 
distance between the macro-asperity contact regions is 
much larger than the linear size ~ Ao/Cm of these re- 
gions). However, owing to the high density of micro- 
asperity contact regions within a macro-asperity contact 
region, it is not possible to neglect the thermal coupling 
(or overlap) between the micro-asperity contact regions 
within each macro-asperity asperity contact region. 

The temperature increase in the macro-asperity con- 
tact regions between a tire and a road surface has re- 
cently been studied using infrared camera. Fig. [5] shows 
a photo of a tire as it rotates out of contact with the road. 
The red-yellow color indicate the "hot" spots arising from 
macro-asperity contacts in the tire-road footprint. 



3 Rubber friction without flash temperature 

The main contribution to rubber friction when a rub- 
ber block is sliding on a rough substrate, i.e., a tire on a 
road surface, is due to the viscoelastic energy dissipation 
in the surface region of the rubber as a result of the pul- 
sating forces acting on the rubber surface from the sub- 
strate asperities, see Fig. El Recently I have developed a 
theory which accurately describes this energy dissipation 
process, and which predicts the velocity dependence of 
the rubber friction coefficient 0, • The results depend 
only on the (complex) viscoelastic modulus E{u>) of the 
rubber, and on the substrate surface roughness power 
spectrum C(q). Neglecting the flash temperature effect, 
the kinetic friction coefficient at velocity v is determined 
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FIG. 3: A rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification (£ = 1) 
it appears as if complete contact occurs in the macro asperity 
contact regions, but when the magnification is increased it is 
observed that in reality only partial contact occurs. 
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FIG. 4: The contact area between an elastic solid with a flat 
surface and a hard randomly rough substrate shown at low 
magnifications (£ = 4), left, and high magnification (£ = 216), 
right. The surface has the fractal dimension Dt = 2.2 and 
qo/qL = 3. Adapted from [Io|. 
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where ao is the mean perpendicular pressure (load di- 
vided by the nominal contact area), and v the Poisson 
ratio, which equals 0.5 for rubber-like materials. The fac- 
tor P{q) = P((qo) = A(()/Aq is the (normalized) area of 
contact when the system is studied at the magnification 

C = q/qo- 




FIG. 5: Infrared photo of tire as it rotates out of contact with 
the road. The red-yellow color indicates the "hot" spots aris- 
ing from macro-asperity contacts in the tire-road footprint. 
Adapted from Ref. [Hill. 



The theory takes into account the substrate surface 
roughness components with wave vectors in the range 
ql < q < Qx, where ql is the smallest relevant wave 
vector of order 2ir/L, where (in the case of a tire) L is 
the lateral size of a tread block, and where gi may have 
different origins (see below). Since qL for a tire tread 
block is smaller than the roll-off wave vector qo of the 
power spectra of most road surfaces (see Fig.|2J), rubber 
friction is very insensitive to the exact value of qL- 

The large wave vector cut off q% may be related to 
road contamination, or may be an intrinsic property of 
the tire rubber. For example, if the road surface is cov- 
ered by small contamination particles (diameter D) then 
qi 2ir/D. In this case, the physical picture is that 
when the tire rubber surface is covered by hard particles 
of linear size D, the rubber will not be able to penetrate 
into surface roughness "cavities" with diameter (or wave- 
length) smaller than D, and such short-range roughness 
will therefore not contribute to the rubber friction. For 
perfectly clean road surfaces we believe instead that the 
cut-off qi is related to the tire rubber properties. Thus, 
the high local (flash) temperatures and high local stresses 
which occur in the tire rubber-road asperity contact re- 
gions (see below), may result in a thin (typically of order 
a few micrometer) surface layer of rubber with modified 
properties (a "dead" layer), which may contribute very 
little to the observed rubber friction (see Sec. 10). Since 
the stresses and temperatures which develop in the as- 
perity contact regions depend somewhat on the type of 
road [via the surface roughness power spectrum C(q)], 
the thickness of this "dead" layer may vary from one 
road surface to another, and some run-in time period 




FIG. 6: The road asperities exert pulsating forces on the slid- 
ing rubber block, leading to energy dissipation in the rubber 
via the internal friction of the rubber. Most of the energy dis- 
sipation occurs in the volume elements bound by the dashed 
lines. The rubber viscoelastic deformations in the large vol- 
ume elements are induced by the large road asperities, while 
the smaller dissipative regions result from the smaller asper- 
ities distributed on top of the large asperities. In general, in 
calculating the rubber friction, the viscoelastic energy dissi- 
pation induced by all the asperity sizes must be included, and 
the local temperature increase (flash temperature) in the rub- 
ber resulting from the energy dissipation must also be taken 
into account in the analysis. 



will be necessary for a new "dead" layer to form when 
a car switches from one road surface to another (such 
"run-in" effects are well known experimentally). 

4 Rubber friction with flash temperature 

The temperature field T(x, t) in the rubber block is 
determined by 



dT , 
— - DV 2 T 
at 



Q(x, t) 
pC v 



(6) 



where Q is the energy production per unit volume and 
unit time as a result of the internal friction in the rub- 
ber. The heat diffusivity D = X/pCy, where p is the 
mass density and A the heat conductivity. For rubber 
we typically have p w 10 3 kg/m 3 , Cy w 10 3 J/kgK and 
A « 0.1 W/mK. This gives D w 10~ 7 m 2 /s. Consider 
now a rubber block with a flat surface sliding on a rough 
hard substrate. Assume that there is an asperity contact 
area with diameter d, see Fig. [7| During sliding at the ve- 
locity v the asperity will generate pulsating forces on the 
rubber surface characterized by the frequency ujq ~ v/d, 
which will result in energy dissipation in a volume ele- 
ment of order d 3 . If the velocity v is high enough, negligi- 
ble heat diffusion will take place during the time period 
d/v and the temperature increase in the rubber in the 
vicinity of the asperity will be AT ps Q/pCy, where Q 
the the amount of frictional heat energy (per unit vol- 
ume) in the volume element d 3 . The assumption of neg- 
ligible heat diffusion requires that the effective contact 



time d/v is smaller than the diffusion time d 2 /D, i.e., 
vd > D. Thus, for example, if v = 0.1 m/s, heat diffu- 
sion will be unimportant if d > D/v ~ 10 -6 m. However, 
after the substrate asperity - rubber contact is "broken" , 
the peak temperature will decrease, and the temperature 
distribution will broaden with increasing time, in a way 
determined by the heat diffusion equation. 

In this section we will develop a theory for the influ- 
ence of the flash temperature on rubber friction when the 
substrate has roughness on very many different length 
scales. We consider first an idealized case where surface 
roughness occurs on a single length scale, and then the 
more complex case of roughness on many different length 
scales. 
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FIG. 7: A rubber block sliding on a rough hard substrate 
with surface roughness on a single length scale. The region 
in the vicinity of a substrate asperity is shown. The heat en- 
ergy production Q(x, t) per unit volume and unit time occurs 
mainly within the volume element surrounded by a dashed 
line. 



We consider a surface with randomness on a single 
length scale Ao = 2-7r/go and with the root-mean-square 
roughness amplitude ho- The surface roughness power 
spectrum C(q) = (hl/2irqo)5(q — go)- The (average) ra- 
dius of curvature of the asperities is R w I/^q/jo)- Fig- 
ure shows the contact between one substrate asperity 
and the rubber. The diameter of the contact area is of 
order d. We assume that no roughness occurs on smaller 
length scale than the size of the asperity in Fig. [7| The 
time it takes to slide the distance d is d/v. During this 
time the heat energy produced in the volume clement 
~ d 3 at the asperity is given by fiF^d, where Fn = od 2 
is the normal asperity load (a is the average perpendicu- 
lar stress in the contact area). Thus the energy produc- 
tion per unit volume Q w pF^d/d 3 = fia. Neglecting 
heat diffusion this will result in the temperature increase 
AT w Q/pCy ~ lAcr/pCy- But we have shown earlier 



ImffQ^o.r) 
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where luq = vqo. From standard contact mechanics the- 
ory one expect a ~ qaho\E(uio,T)\. Thus the friction 
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lmE(u ,T) 
\E(w ,T)\ 
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and the temperature increase 



or 



AT a (q h ) 



T « T + (q h ) 



2 ImE(uJo,T) 



pC v 
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where To is the background temperature. Note that the 
complex elastic modulus E(ui,T) depends on the (lo- 
cal) temperature T. For "simple" (unfilled) rubber the 
Williams- Landel- Ferry equation ( WLF) [l3| equation can 
be used to (approximately) describe the temperature de- 
pendence of E(u>,T): 



E(lo,T) = E{ua,T / &t , To) 



where 



T — T 

log 10 a T » —8.86 
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51.5 + T-T„ 
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For any given viscoelastic modulus E(lj,Tq), equations 
(7)-(9) form a complete set of equations from which the 
temperature T and the friction coefficient p can be ob- 
tained by, e.g., iteration. Here I will only discuss quali- 
tatively how the temperature increase AT influence the 
rubber friction. 

In order to understand the qualitative influence of the 
flash temperature on rubber friction it necessary to know 
the general structure of the viscoelastic modulus E{u) 
of rubber-like materials. In Fig. [5] we show the real 
E\ = RcE and the imaginary part E 2 = ImE of E(ui) 
and also the loss tangent E 2 /Ei. At "low" frequencies 
the material is in the "rubbery" region where KeE(uj) is 
relative small and approximately constant. At very high 
frequencies the material is elastically very stiff (brittle- 
like). In this "glassy" region ReE is again nearly constant 
but much larger (typically by 3 to 4 orders of magnitude) 
than in the rubbery region. In the intermediate frequency 
range (the "transition" region) the loss tangent is very 
large and it is mainly this region which determines, e.g., 
the friction when a tire is sliding on a road surface. 

The influence of the flash temperature on rubber fric- 
tion differs depending on if the perturbing frequency 
luq = vqo is smaller or larger than the frequency uj\ where 
Ym.E(bj,Ta)/\E(u>,To)\ is maximal. Since an increase in 
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FIG. 8: (a) The viscoelastic modulus E(uj) = Ei + iEi of a 
typical rubber-like material, and (b) the loss tangent E2/E1. 
(Schematic.) 



the temperature (To — > T > To) shifts the viscoelastic 
spectrum to higher frequencies (see the WLF function), 
if luq < u>i, the flash-temperature will decrease the fric- 
tion. On the other hand, if luq > w\ the opposite effect 
occur. However, in most applications, e.g., in tire appli- 
cations, the perturbing frequencies are (almost) always 
below u>i and the friction will decrease when the flash 
temperature is taken into account. This has extremely 
important practical consequences, as we will be discussed 
later. 

We consider now the role of the flash temperature for 
the general (and more complex) case where roughness 
occurs on many different length scales. We consider first 
stationary sliding and then non-stationary sliding. 

4.1 Stationary sliding 

In this section we will develop a general expression for 
the friction acting on a rubber block sliding at a constant 
velocity on a randomly rough substrate. We will take into 
account the effect of the flash temperature. We will in- 
cluding the thermal overlap between the heat produced 
in the micro-asperity contact area inside every macro- 
asperity contact area. That is, the temperature rise at 
one micro-asperity contact area will produce a subse- 
quent temperature rise at a neighboring micro contact. 
We will make use of a "mean field" type of approximation 
where we laterally smear out the heat sources associated 
with the micro-asperity contacts within a macro-asperity 
contact, which should be an excellent approximation in 



*: : : 
* 

/ 

* : . * . * : * > .*.*.* 




. . % 


V 


ri 




\ - ': 
: :\: : 

: i : V 

• . • . \ \ ; ; 

• • • \ 




. ; .........v-v. v • . 1 : ■ ; . ; 

r / 

/ ■ 1 ~ V 

/ V T 2 ' 

«-v' V A '>' 




\ 

" "■ \ ■" " " 
\ 

^ 1 

; 




. . .. ■ 


di 


V 1 

d 2| 






w 











FIG. 9: A rubber block sliding on an asperity with shorter 
wavelength asperities on top of it. The temperature increases 
T < T < T 2 . 



most cases. 

Let us now discuss the flash-temperature when a rub- 
ber block is sliding on a hard substrate with surface 
roughness on many different length scales. The basic 
problem is illustrated schematically in Fig. [5] in the case 
of only two length scales, where a "large" asperity is cov- 
ered by "small" asperities. In response to the large as- 
perity there will be a heating of the rubber in the big 
volume element of linear size d\ . If To denote the back- 
ground temperature, then the flash temperature in the 
large volume element will be higher, say T\ = T + ATi . 
The flash temperature in a small-asperity volume ele- 
ment (linear size d,2, see Fig. will be even higher, say 
T2 = Ti + AT2. For surfaces with surface roughness on 
many length scales the temperature will increase as we 
go to smaller and smaller asperity contact regions. We 
now study the temperature distribution quantitatively. 

The average shear stress which acts in the macro 
asperity-contact area is 



An 



(Tf 



(10) 



where at is the nominal frictional shear stress, and A(£ m ) 
the macro asperity contact area observed at the magni- 
fication C = C, u = Ao/A m = q m /qo which usually is of 
order unity (see Appendix A). Using (3) and (10) and 
that P(q) = A(()/A gives 
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P(q) E{qv cos <f>,T q ) 



Im 



1 



(11) 

where T q is an effective temperature in the volume of the 
rubber involved in the calculation of the contribution to 
rubber friction from surface roughness on the length scale 
A = 27r/g (see below). Thus the energy "dissipation" per 



unit time and unit area in a macro contact area 
T « f 2 2 , n( x p (g) j cos ^,T g ) 

(12) 

The energy production per unit volume, from asperities 
on the length scale A = 27r/q, decay into the solid as 
r~j e~ 2qz so that the energy production per unit volume 
and unit time in a contact area is obtained by introducing 
in the integrand in (12) the factor 
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Using this equation the energy production per unit vol- 
ume and unit time in the rubber in a macro-asperity 
contact area is given by: 
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where we have assumed that the energy production start 
at t = 0. We write the step function 8(t) as 

6(t)=~ [ duj — ^e"-* 
2m J oj + %e 

where e is a positive infinitesimal number. We can also 
write 



Thus 
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where 
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(16) 

To solve the heat diffusion equation (6) we need a bound- 
ary condition on the surface z — 0. We will assume that 
dT/dz = for z — 0, i.e., we neglect heat transfer from 
the rubber to the substrate. This is an excellent approx- 
imation at high enough sliding velocity. At low sliding 
velocity it will give rise to an overestimation of the flash 
temperature, but in this case the flash temperature effect 
is anyhow not very important. The boundary condition 
dT/dz(z — 0) = is equivalent to solving the heat diffu- 
sion equation for an extended solid (—00 < z < 00) with 



a symmetric heat source obtained by replacing exp(— 2qz) 
in (14) by exp(— 2g|z|). Thus we get 



T{ Zl t) = T + j 
Jo 



dq f(q) 



X-7 TTT / <LC(1J; 

(2tt 



1 



1 



oj + iDk 2 oj + ie k 2 + 4q 2 



Performing the ^-integration gives 

T(z,t) = T + dqf(q) 
Jo 



— / '/A- — ( 1 - c- DkH \ ,„ 4g , „ e ikz 



(17) 



2tt 



Dk 2 



k 2 + Aq 2 



(18) 



The optimum (spatially uniform) temperature to be used 
when calculating the contribution to the friction from 
surface roughness on the length scale A = 2n/q can be 
obtained using 



T a = 



f °°dz T(z,t)e- 2 ^ 



f °°dze-M- 
Using (18) this gives 

T q = T Q + / dq' g(q, q')f(q') 
Jo 

where 

1 r°° 1 / 

9M = -J Q dk Dk^ 1 



(19) 



(20) 



— e 



Dk 2 t 



Aq' 



Aq 2 



k 2 + Aq' 2 k 2 + Aq 2 



(21) 



where to f=a R/v is roughly half the time a rubber patch 
is in contact with the macro asperity. The radius R of a 
macro-asperity contact region is estimated in Appendix 
A. Note that the complex elastic modulus E(ui,T) de- 
pends on the (local) temperature T. In general one may 
write E(u>,T) — bTE(u)aT,To), where or and &t depend 
on the temperature T but with ax — bT = 1. For 
unfilled rubber br ~ 1 and ax is (approximately) given 
by the WLF equation. The functions clt and 6t are to- 
day measured routinely using standard rheological equip- 
ment. 

When T q has been obtained from (20) and (21), the 
friction coefficient can be calculated using the equations 
derived in Ref. |9j: 



1 



in 



f, a - / dq q* C(q)P(q) 
L Jqo 



COS 



Im 



E(qv cos </>, Tq) 
(l-^>o 



P(q) = 1 (°° dx ^exp [-x 2 G(q)} 
n Jo x 



where 



G(q) = l f dqq 3 C(q) f d<P 

8 Jqo Jo 



E(qv cos </>, T q ) 



(1 - ^>o 



4.2 Non-stationary sliding 

Here we develop a completely general theory of non- 
stationary sliding. Consider again the heat diffusion 
equation 



(ii-|x-r(t)|) (22) 



where v(t) — xx(t) is the position vector of the bottom 
surface of the rubber block, and where R is the radius of 
a macro-asperity contact region which we assume to be 
circular. Using (14) and (16) gives 



pC v 



= dq f(q, 
Jo 



t) e 



-2qz 



(23) 



Let us introduce a new coordinate system, moving with 
the bottom surface of the rubber block: 

x' = x - r(t) 

Substituting this in (22), and replacing x' — > x for sim- 
plicity, gives 

dT .,, dT d 2 T Q{z,t) n , „ , N 

We now introduce the Fourier transform 

f(p, k, t) = ^ J d 2 xdz T(x, z, t)e-^ +kz ^ (25) 



T(x,z,t) = J d 2 pdk f(p,fc,t)e i(p ' : 



:+fcz) 



(26) 



Using (25) we get from (24): 



^-zp-r(t)T + Dk 2 T 



Q(k,t) 1 f R 



drrJ (kr) 



(27) 



pC v 2ir J 

where Jo (a;) is the zero order Bessel function and where 
Q(k, *) = ^ Jdz Q(z, t)e- lkz . (28) 
Equation (27) is easy to integrate to get 



f(p,k,t) = [ dt' 
Jo 



-Dk 2 (t-t')+ip :c [x(i)-x(t')\ 



Q(k,t') 1 

pC V 27T Jo 



drrJo(kr) 



(29) 



The laterally averaged temperature in the contact area 
is 



T(M) 



1 



kR 2 J\ x \<r 



d 2 x T(x,fc,i) 



\x\<R 



d 2 x J d 2 p e ipx T(p, k, t) 



= jpjd 2 p J drrJ (kr)f(p,k,t) (30) 
Substituting (29) in (30) gives after some simplifications 

f(k,t)= f dt 1 T(t,t')e- Dk2 ^^^- (31) 
Jo P^v 



where 



2 f°° 

W) = / dppj (p[x(t) - x(t')]) 



(L 



x / drr.Jo(pr) 



Now since 



/ drrJ {pr) = —Ji(pR), 
Jo P 

we get with pR = y 

f-OO 

T(t,t') = 2 / dyy- 1 J Q (2wy) (Ji(y)) 2 = h(w) (32) 
Jo 

where 

w = w(t, t') = [x{t) - x(t')]/2R (33) 
Performing the integral in (32) gives 



h(w) = 1 - -w (l - w 2 ) 1/2 - -arcsin(w) (34) 

7T 7T 



for w < 1 while h = for w > 1. Using (23) we get 
Q(M) 1 f° j_ ^ 4g 



sj£ *'<«■•>**£? (35) 



pCy 2tt,/o ™" J K ' 1, " / k 2 +4q 2 
Next, using that 



/>0O 

fg(t) = 2q / f(z,t)e 
Jo 



Iqz 



•■I 

Jo 



2 \ dk 



4q 2 
k 2 + Aq- 



:f(k,t) 



(36) 



we get from (31), (32), (35) and (36) 



f q (t) =T 



f dt' T(t,t') f 
o Jo 



dq'f(q',t') 



In the limit D — > oo we get as expected T — To- 

Finally, let us consider stationary sliding so that x(t) = 
v t. In this case w = (v /2R)(t — t') and f(q,t) = f(q) 
is independent of time. Thus (39) takes the form 



T Jo 



x- / dk 



fc 2 + V k 2 +4q' 2 



(37) 



Let us summarize the basic equations derived above. 
The friction coefficient 



\ C 'dqq 3 C{q)P{q,t) 



f 27T , , , T E(qv(t) cos 4>,T q (t)) 

x / dej) cos 4> Im w ) ; gv yy (38) 

Jo (1 - ^ jo-o 

where ctq — F/v(i)/Ao, where Ao is the nominal contact 
area. In this equation enters the flash temperature at 
time t: 



f-t f-OO 

T q (t)=T Q + / dt'T(t,t') / dq'f(q',t' 
Jo Jo 



7T J Q k 2 + 4q 2 k 2 +4q' 2 



(39) 



where T(t,t') = h(w(t,t')) with w{t,t') = \x{t) - 
x(t')]/2R depend on the history of the sliding motion. 
The function 



f{q,t) = ^q A C{q) 



pC v P(q m ,t) 



I 



x / d<f> coscj) Im 



E(qv(t) cos 4>,T q (t)) 



\-v 2 



(40) 



where v = x(t) depends on time. The function P(q,t) 
(which also depends on time) is given by the standard 
formula: 



, . 2 f°° smx 

P(q,t) = - dx — 

7T Jo X 



where 



exp [— x 2 G(q, t)] = erf 



1 



2^G 
(41) 



G{q,t) = \ f q dqq 3 C(q) f* deb 

° Jqo Jo 



E(qv(t) cos 4>,T q (t)) 



(1 - ^>o 



(42) 



4.2.1 Limiting cases and physical interpretation 

of h(w) 

When the heat diffusivity D = equation (39) takes 
the form 



ft fOC 

f q (t)=T + / dt'T(t,t') / dq'f(q',t' 
Jo Jo 



>iT7 (43) 



T q (t)*T + dq'G(q,q')f(q') (44) 



where 



Aq 2 



tt J k 2 + 4q 2 k 2 +4q' 



T M(k) (45) 



where 



M(k) 



dt' h(w)e 



(46) 



>t-2R/v 

Changing the integration variable from t' to w gives 

M (k) = — du> ^( w ) e -( 2 «/^) Dfc2t " (47) 
«o Jo 

The result (45) is the same as (21) but with a slightly 
different (and more accurate) function M(k). Thus, in 
(21) occur instead of M(k) the factor 

^2 (l - e- Dfe2 *°) = t dw e - toDk2w 

Since to is of order R/vq and h(w) of order 1, these two 
expressions for the steady rubber friction are of similar 
magnitude. 




FIG. 10: The common region (dotted area) between the two 
circles has the area nh{w). 

Here I present an alternative derivation of the factor 
h(w) = g(t,t') which gives a direct physical interpreta- 
tion of this factor. In the theory above this function 
emerged purely mathematically, as a result of an integral 
of a product of three Bessel functions. 

Let us assume that the heat diffusivity D — and 
consider the equation 

^=K(z,t)0(R-\x-T(t)\) (48) 

where K = Q/pCy and where 9(x) is the step function. 
Here R is the radius of the macro contact area and r(t) 



the position vector of the contact area. Integrating (48) 
gives 

T(x,z,t)=T + f dt' 6{R- |x-r(t') \)K(z,t') (49) 
Jo 

Averaging the temperature (49) over the contact area 
gives 

f(z,t)=T + -^ [ d 2 x 

J\ x -r(t)\<R 

x / dt' 0(R-\x-r(t')\)K(z,t') 
Jo 

i r* 

= T ° + -m dt ' K (*J) 

x[ d 2 x 6{R- |x + r(t) -r(t')|) (50) 

J\x\<R 

If we change integration variable x' = x/i? and denote 
x' with x for simplicity, we get from (50) 

T = T + - / dt' K(z,t') [ d 2 x(9(l-|x+2w|) (51) 
T Jo '|x|<l 

where 

w = [r(t) - r(t')]/2R 
Note that the integral 

h(w) = - I d 2 x 6(1- |x + 2w|) (52) 

n J\x\<l 

is the common area (divided by n) between two circular 
regions (with unit radius) with the origins separated by 
the distance 2w, see Fig. It is easy to calculate 

2 1/2 2 
h(w) — 1 w (l — w ) arcsin(w) (53) 

7T 7T 

for w < 1 while /i(w) = for io > 1. This formula agrees 
with (34). Note that w = w(t,t') and = r(t, t') so 
that (51)-(53) gives 

1 I* 

T = T + - / dt' r(i, t')K(z, t') (54) 

Finally, using (14), (16) and (19) in (54) gives 

f q (t)=T + [ dt' T(t,t') [°° dq'f(q',t')^ 1 
Jo Jo 1 + 1 

which is the same as (43). Thus, we conclude that the 
factor h(w) — T(t, t') in (37), is equal to the overlap area 
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FIG. 11: The friction coefficient (a), the flash temperature 

(b) , and the logarithm of the (normalized) area of contact 

(c) , as a function of the logarithm of the sliding velocity. For 
a rubber block (tire tread rubber) sliding on an asphalt road 
surface. The background temperature To = 60 °C and ( mic( = 
400. In (a) we show the kinetic friction coefficient both with 
and without the flash temperature. In (b) and (c) we show 
results both for f = 1 and at the highest magnification £ = 
400. 



(divided by w) between the contact area at time t, and 
the contact area at an earlier time t', see Fig. 

5 Sliding dynamics: numerical results 

In this section we will present numerical results for the 
friction force acting on a rubber (tire tread) block, slid- 
ing on an asphalt road surface. We will study the friction 
as a function of the velocity of the bottom surface of the 
block. In real experiments it is, of course, not possible to 
specifies the motion of the bottom surface directly (unless 
steady sliding occur) , but usually one specify the motion 
of the top surface of the rubber block (or some other dis- 



tant part of the system). We will consider this case later. 
In this section we assume the rubber background temper- 
ature T = 60 °C, the substrate is an asphalt road surface 
(with the power spectrum given in Fig. |2J) and include 
the substrate roughness down to the wave vector cut-off 
Qi = Cmax<7o with Cmax = 400 and qo — 1500 m _1 so that 
qi = 6 x 10 5 or Ai = 2n/qi w 10 /im. All calculations is 
for the nominal (squeezing) pressure <tq = 0.4 MPa. 
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FIG. 12: The friction coefficient (a), the flash temperature 

(b) , and the logarithm of the (normalized) area of contact 

(c) , as a function of the logarithm of the magnification. For 
a rubber block (tire tread rubber) sliding on an asphalt road 
surface. The background temperature To = 60 °C and the 
sliding velocity v — 1 m/s. In (a) and (c) we show results 
both with and without the flash temperature. 



5.1 Stationary motion 

Fig. II II shows the friction coefficient (a), the flash tem- 
perature (b), and the logarithm of the (normalized) area 
of contact (c), as a function of the logarithm of the slid- 
ing velocity. In (a) we show the kinetic friction coeffi- 



cient both with and without the flash temperature. In 
(b) and (c) we show results both for £ = 1 and at the 
highest magnification £ = 400. Note that without the 
flash temperature the rubber friction coefficient increases 
monotonically up to the highest studied sliding velocity 
v w 100 m/s. When the flash temperature is included in 
the study the friction is maximal for v ps 1 cm/s. The 
decrease of the friction observed when the flash tempera- 
ture is included in the analysis is easy to understand (see 
also Sec. 4): when the rubber heats up the viscoelastic 
modulus E{u) shift towards higher frequencies and the 
rubber becomes more elastic (less viscous) resulting in 
less energy dissipation. 

Fig- Hir b) shows that the flash temperature at the 
magnification £ = 400 (i.e., the temperature about 10 /im 
below the rubber surface) is « 140 °C, i.e., about 80 °C 
above the rubber background temperature. On the other 
hand the temperature increase a few mm below the sur- 
face (corresponding to the magnification £ = 1) is just 
a few degree. We note that during steady sliding the 
temperature in the whole rubber block will increase con- 
tinuously with increase time, but this effect is not in- 
cluded in the discussion above, unless one allowed for 
the background temperature T to increase with increas- 
ing time. The calculation of the time dependence of 
To require the knowledge of the temperature of the sur- 
rounding medium, and the heat transfer coefficient to the 
surrounding. This topic is important in many practical 
applications, but does not interest us here. 

Fig. Illf c) shows that the area of (apparent) contact 
at low magnification (£ = 1) decreases continuously with 
increasing sliding velocity. This result is expected be- 
cause the frequencies oj of the pulsating deformations of 
the rubber increases with increasing sliding velocity v 
(lo ~ i>), and rubber becomes elastically stiffer at higher 
frequencies, resulting in an (apparent) smaller contact 
area at high sliding velocity. However, at the highest 
magnification £ = 400, fig. Illf c) shows that the area 
of contact increases with increasing sliding velocity for 
v > 0.1 m/s. This is caused by the increase in the 
temperature in the surface region of the rubber, which 
shift the viscoelastic modulus towards higher frequencies 
faster than the linear increase to ~ v in the perturbing 
frequencies arising from the increase in the sliding veloc- 
ity. 

Fig. 1121 shows (a) the friction coefficient, (b) the flash 
temperature, and (c) the logarithm of the (normalized) 
area of contact, as a function of the logarithm of the 
magnification, for the sliding velocity v = 1 m/s. In 
(a) and (c) we show results both with and without the 
flash temperature. Fig. I12f a) shows that when the flash 
temperature is included, the surface roughness on every 
decade in length scale contributes roughly with an equal 
amount to the rubber sliding friction. Fig. I12f c^l shows 
that at low magnification the contact area is the same 
with and without the flash temperature included in the 



analysis, while for large magnification the area of contact 
increases when the flash temperature is included, as ex- 
pected from the elastic softening of the rubber at higher 
temperatures (see above). 
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FIG. 13: Non-uniform sliding motion, (a) The velocity of 
the bottom surface of the block increases linearly with time 
from zero to 1 m/s and then decreases back to zero. In case 
1 the time for the whole velocity cycle is 0.002 s and for case 
2, 0.004 s. (b) The friction coefficient exhibit hysteresis as a 
function of the sliding velocity, and is smaller during the retar- 
dation time period than during the acceleration time period. 
This behavior reflect the build up of the flash temperature 
during the sliding cycle. 

5.2 Non-stationary motion 

It is usually assumed in simple treatments of friction, 
that the friction coefficient is a function only of the 
instantaneous slip velocity Vh(t), i.e., = /i(vb(£))j 

where Vb(t) is the velocity of the bottom surface (in con- 
tact with the substrate) of the sliding block. However, 
this approximation fails badly for rubber friction. The 
reason is the strong dependence of rubber friction on the 
temperature distribution in the rubber block. Since the 
temperature distribution T(x, t) at time t depends on 
the sliding motion for all earlier times t' < t, it follows 
immediately that for rubber-like materials the fiction co- 
efficient will depend on the sliding history, i.e., it will be 
a functional of Vb(t): [i — /j,(vb(t'),t' < t). In this section 
I present numerical results which illustrate this fact. 
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FIG. 14: The flash temperature (a) and the frictional shear 
stress (b) as a function of time as the bottom surface of the 
rubber block perform the motion indicated in Fig. HSf a). and 
with u b = for t > 0.002 s and 0.004 s for case 1 and 2, 
respectively. 



Let us assume the velocity of the bottom surface of 
the block first increases linearly with time from zero to 

1 m/s and then decreases back to zero, as indicated in 
Fig. I13f a). In case 1 the time for the whole cycle is 
0.002 s and for case 2, 0.004 s. In Fig. Hlb) we show 
the resulting friction coefficient. Note that fi exhibits 
hysteresis as a function of the sliding velocity, and that 
// is smaller during the retardation than during the ac- 
celeration time period. This behavior just reflects the fi- 
nite slip-distance necessary for building up the flash tem- 
perature, and the finite time involved in heat diffusion, 
and shows that (i does not just depend on the instan- 
taneous sliding velocity but on the whole sliding history 
(memory-effects). Note that the slip-distance is given by 
n^/o, where u ma x = 1 m/s is the maximum of the ve- 
locity and the acceleration a — 1000 and 500 m/s 2 for 
case 1 and 2 respectively, giving the slip distances 1 and 

2 mm, respectively. These values are both smaller than 
the diameter of the macro-asperity contact regions which 
is 2R « 3.3 mm. The larger slip distance for case 2 im- 
plies that more energy is deposited (in the rubber) at the 
rubber-substrate contact regions than for case 1, which 
will result in a higher flash temperature in case 2 than 
in case 1. This explains why the friction coefficient is 
smaller for case 2 (see Fig. I13f b'lb 
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FIG. 15: (a) the friction coefficient, (b) the flash-temperature 
at the highest magnification and (c) the (relative) contact 
area at the highest magnification, as a function of time. The 
bottom surface of the block moves with the constant velocity 
v h = 2 m/s for < t < 0.002 s, and at t = 0.002 s it is 
abruptly reduced to 10~ 5 m/s and is kept at this velocity for 
0.0005 s (case 1), 0.001 s (case 2), 0.002 s (case 3), 0.003 s 
(case 4), and then returned to Vb = 2 m/s. 



The temperature buildup during the sliding motion is 
illustrated in Fig. This figure also shows that after 
the sliding motion has stopped, i.e., for t > 0.002 s and 
t > 0.004 s for case 1 and 2, respectively, the temperature 
at the surface decreases monotonically because of heat 
diffusion. Note that as a function of the slip velocity 
«b (not shown), the flash temperature in Fig. I14f a % ) is 
always higher for case 2 as compared to case 1. Fig. 
I14f b) shows the frictional shear stress <J[(t) — <7o/x(i) as 
a function of time. 

Assume now that the bottom surface of the block 
moves with the constant velocity = 2 m/s for < 
t < 0.002 s, and at t = 0.002 s it is abruptly reduced to 



10 5 m/s, and is then kept at this value for 0.0005 s (case 
1), 0.001 s (case 2), 0.002 s (case 3), 0.003 s (case 4), 
and then returned to — 2 m/s. For this sliding history 
we show in Fig. ^| (a) the friction coefficient, (b) the 
flash-temperature at the highest magnification and (c) 
the (relative) contact area at the highest magnification, 
as a function of time. Note that the longer the system is 
kept in the low-velocity state (where v = 10 -5 m/s), the 
higher the "start-up" peak in the friction will be when 
the velocity is switched back to 2 m/s. The reason for 
this is that, because of heat diffusion, in the low- velocity 
state the temperature in the rubber will decrease con- 
tinuously with increasing time, see Fig. I15f b). and the 
rubber friction increases when the rubber temperature 
decreases. Fig. E3 C ) illustrate that, as expected, the 
contact area between the rubber and the substrate in- 
creases when the rubber heats up and when the rubber 
sliding velocity decreases. 




FIG. 16: A rubber block (thickness d) pulled on a rough hard 
substrate. The upper surface is clamped and moved with a 
prescribed velocity v(t) = x t . In the mathematical descrip- 
tion of the sliding motion, the block is discretized into N 
layers of thickness Ad = d/N. The coordinates xi(t) = Xb{t), 
X2(t), x^(t),xn+i = x t (t) depend on time. The continuum 
limit is obtained as N — » oo. 



6 Rubber block dynamics 

We consider the simplest but most fundamental case of 
a rubber block (mass M) squeezed against a substrate, 
and with the upper surface clamped and in prescribed 
lateral motion as indicated in Fig. 

The shear stress that acts on the lower surface of the 
rubber block is treated as a spatially uniform stress <Tf(t). 
In the most general case, where inertia effects are impor- 
tant, we discretize the block (thickness d) into N layers 
of thickness Ad = d/N and mass m = M/N. The coor- 
dinates xi(t) — Xb(t), Xi{t), xu(t), xa+i(t) = x t (t) 
depend on time. The continuum limit is obtained as 
N — > oo. The non-stationary sliding of the rubber block 
[see Fig. I16j is determined by the equations of motion 



mx b (t) = Fx(t) - F [x b (t)}, 



mx 2 (t) = F 2 {t) - Fx(t), 



mx 3 (t) = F 3 (t)-F 2 (t), 



or, if we use (57), 



mx N (t) = F N (t) - F N -!(t), 



(55) 



where Fq is the friction force acting on the bottom surface 
of the block, and Fi the shear force acting on layer i from 
the layer i + 1 above. Fn is the force acting on layer N 
from the drive (the black slab on top of rubber block in 

Fig. inj. 

Note that the friction force Fo[xi,(t)] is a functional of 
Xf,(t) since it depends on the whole history {xb(t');t' < 
t}. The shear force Fi acting on layer i from the layer 
above, can be derived from the equation relating the 
shear stress o~{u>) to the shear strain e(u>) via 



ct(oj) = G(aj)e(uj), 



(56) 



where G{lo) is the complex frequency dependent shear 
modulus. In our case (i — 1, N) 



s(t) = [x i+1 (t)-Xi(t)]/Ad, 



(57) 



where Ad is the thickness of the rubber block. If we 
multiply (56) with the area Aq = L x L y of the tread block 
we get the shear force 



Fi(u) = A,G(w)eH, 



or in time-space: 



Fi(t) = / du Fi(u)e 



-ILOt 



du A()G(uj)e(uj)e~ 



(58) 



Substituting 



e(o;) = ^Jdt' e(t')e^' 
in Eq. (58) gives 

/oo 
dt' A G(t-t')e(t'), 
-oo 

where 

1 f°° 

G(t) = — dcu G(u)e 

27r J-oo 



-iujt 



(59) 



(60) 



If we assume that e(t) = for t < and if we use that 
G(t) must vanish for t < because of causality (see Sec. 
7), we get 



Fi(t)= / dt'A G(t-t')e(t'), (61) 
Jo 



Fi(t)= dt' K(t-t')[x i+1 (t')- Xi (t% (62) 
Jo 

and in particular 



F N (t)= / dt> K(t-t')[x t (t')-x N (t% (63) 
Jo 

where the memory spring 

K(t) = A G(t)/Ad. (64) 

The equations above can also be applied when the 
properties of the rubber block depends on the vertical co- 
ordinate z if one replace the mass m — > m, and memory 
spring K — > JQ. If p = p(z) is the rubber mass density, 
then rrii = AoAdp(zi) and similarly Ki = AoG(zi, t)/ Ad. 
Since the rubber friction Fq depends extremely strongly 
on the sliding velocity for very small sliding velocities, 
the system of differential equations given above is of 
the stiff nature, which requires special care in the time- 
integration in order to avoid numerical instabilities. 

7 Viscoelastic modulus 

Many experimental techniques can be used to obtain 
the complex, frequency dependent, shear modulus G{u) 
[or Young modulus E(u)]. This quantity enters directly 
in the calculation of the rubber frictional shear stress 
Of. But in the rubber block dynamics theory described 
in Sec. 6 we also need the time-dependent real shear 
modulus 



G(t) 



1 

2^ 



du> G(uo)e 



(65) 



Since G(w) usually is known only in a finite frequency 
interval, luq < u> < ui\, it is not easy to calculate G(t) 
directly from (65) using, e.g., the Fast Fourier Transform 
method. Another, even more serious problem is that in 
some cases G(t) is not a "perfect" linear response func- 
tion (see below), which gives rise to serious problems if 
one tries to calculate G(t) directly from (65), e.g., one 
finds that G(t) is complex rather than real. 

Let us first assume that G(cu) is the viscoelastic mod- 
ulus of a linear viscoelastic solid. Now, the shear stress 
o~(t) in a solid at time t can only depend on the deforma- 
tions (or strain) the solid has undergone at earlier times, 
i.e., it cannot depend on the future strain [causality). 
Thus we can write 

a(t) = f dt' G(t-t')e(t'). 



The Fourier transform of this equation gives 
a~(u>) = G(u))e(u>) 
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FIG. 17: The real and the imaginary part of the viscoelastic 
modulus G(u>) for Styrene Butadiene copolymer (a) without 
filler, and (b) with 30% carbon black filler. The red and 
green lines are the original data and the fit-data obtained by 
fitting the experimental G(uj) to a sum over relaxation times 
as described by (68). For the temperature T — 60 °C. 



where 

poo 

G(u) = dt G{t)e iuJt 
Jo 

Since Re(iajf) < for t > and lm(u>) > it follows that 
G(uj) is an analytical function of lo in the upper half of 
the complex frequency plane. Thus all the singularities 
(poles and branch cuts) of G(u>) must occur in the lower 
part of the complex aj-planc and we may write 



G(w) = G 



dr 



H(r) 
1 — ilot 



(66) 



where the spectral density H{t) is a real positive function 
of the relaxation time r. This representation of G(u>) has 
the correct analytical properties (all the singularities oc- 
cur in the lower part of the complex w-plane, as required 
by causality). The idea is now to determine G x and 
H{(J) so that the difference AG(w) = G mcas (ui) — G(ui) 
between the measured G mea s shear modulus and the an- 
alytical expression G{uS) given by (66) becomes as small 
as possible. Thus, we minimize the quantity 



V 



1 



w Jujq |LJmcas(^)| 



(67) 



Note that 

AG = G m eas(w) — Goo + 



H(r) 



dr- . 

1 — ILOT 



In practice, G mca s is only known at a set of discrete fre- 
quencies u> n , n = 1, 2, N. Thus, (67) must be replaced 
with 



V 



1 N 

-Y 



|AGK)| 

N^\G (0J n )\ 



Furthermore, in numerical calculations it is only possible 
to include a finite number of relaxation times in the 
spectral representation of G(w) and write 



GH « Goo - J2 T 



lUJT k 



(68) 



where 



If the relaxation times are distributed exponentially (i.e., 
T k = T o ex P(%)j where tjk are uniformly distributed), it 
is usually enough to include ^10—15 relaxation times, 
or one relaxation time for each decade in frequency. If 
the measured G mea s(w), luq < lo < u>i, correspond to 
a true linear response function then one can find Goo 
and H(lo) so that V vanish. Since in many cases G(w) 
is not a true linear response function (see below) it is 
impossible to choose the real function H(u>) so that V 
vanish. However, if we can find H(oj) and Goo so that 
V « 1 we have a very good representation of G(u>) of 
the form (66). Using (66) it is easy to integrate (65) to 
obtain 



G{t) = Goodit) - 6{t) 



dr m e -*/r 



or, for a finite set of relaxation times, 



G{t) = Goo5{t) 



k Tk 



The minimization of V with respect to Goo and Hk 
can be performed by using a random number generator: 
First some (arbitrary) set of parameters Goo and Hk are 
chosen and V calculated. Next, one replace Goo and Hk 
with Goo(l + £) and Hk(l + ^k), where £ and £k are small 
random numbers. If the new V , calculated using these 
parameter values, is smaller than the original one, then 
the new parameters are accepted as improved parame- 
ters; otherwise they are rejected and the old parameters 
kept. If this procedure is repeated (iterated) many times 
the effective "potential" V will converge towards a min- 
ima, and the parameters Goo and Hk at the minima of 
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FIG. 18: The spectral density as a function of the re- 
laxation time Tfc as obtained by fitting the measured G(w) 
to the spectral representation (68). The squares and circles 
correspond to the unfilled and filled SB copolymer. For the 
temperature T = 60 °C. 

V are the best possible ones, which we use in calculating 
G(t) as described above. 

As an example, in Fig. I17f a) we show the real and 
the imaginary part of the viscoelastic modulus G(o>) for 
Styrene Butadiene copolymer (a) without filler, and (b) 
with 30% carbon black filler. The red and green lines 
are the experimental data and the fit-data, respectively, 
where the fit-data was obtained by approximating the 
experimental G(uj) by the sum (68) as described above. 
Note that for the unfilled rubber both the real and imag- 
inary part of G(uj) are very well fitted by suitable choice 
of a single real (positive) function H(t). This is possible 
only because ReG(w) and ImG(w) are not independent 
but related via a Kramers-Kronig relation (see below). 
For the filled rubber the fit of ReG is equally good, but 
there is a small deviation for ImG. We have observed 
the same for other rubber compounds, and sometimes 
the deviation AG between measured and fitted ImG is 
larger than in Fig. I17f b). but the calculated G{t) is 
still accurate enough for our applications. The devi- 
ation AG reflects the fact that for filled rubber G(uj) 
is not a "perfect" linear response functions, but exhibit 
non-linearly. One source of non-linearity is the so called 
Payne effect, due to the strain- induced break-up of the 
filler networkdHIU. 

Fig. ^] shows the spectral density Hk as a func- 
tion of the relaxation time rfc as obtained by fitting 
the measured G(cu) to the spectral representation (68). 
The squares and circles correspond to the unfilled and 
filled SB copolymer. Note that the main difference be- 
tween the two cases occur for the longest relaxation times 
r > 10~ 8 s, where the filled rubber has a much higher 
spectral density H(t). This implies that for the filled 
rubber, in contrast to unfilled rubber, there is a high spa- 
tial density of (relatively) high-energy-barrier rearrange- 
ment processes. This results from polymer molecules 



bound to the surfaces of the carbon filler particles; the 
bound layers, which may be ~ 1 nm thick, are in a glassy- 
like state and exhibit larger energy barriers for (thermally 
activated) rearrangements, as compared to the polymer 
segments further away from the filler particles. 

We emphasize that only when G(ui) corresponds to a 
linear response function causality will result in a well- 
defined analytical structure for G(cu). In this case the 
real part Gi = ReG and imaginary part Gi = ImG of 
G(iS) are related via Kramers-Kronig relations: 

G 1 (u)-G ao = -P f°° (69) 
G 2 M = -ip r dJ GlK) ' Goc , (70) 

where G^ = G(oo) is the (real) high frequency shear 
modulus. For filled rubbers, (69) and (70) only hold ap- 
proximately. In particular, due to the Payne effect (see 
above) the rubber behaves in a non- linear way i.e., G{u>) 
is not a "perfect" linear response function, and eqs. (69) 
and (70) will only hold approximately. 

8 Numerical results: rubber block dynamics 

In this section we present results for the sliding dy- 
namics of a rubber block with thickness 0.5 cm, squeezed 
against an asphalt road with the nominal pressure cto = 
0.4 MPa. 

Assume first that the upper surface of the rubber 
block is clamped and undergoes uniform acceleration 
a = 500 m/s 2 for 0.004 s. Fig. H3 shows (a) the ve- 
locities v t and «b of the top and bottom surface of the 
rubber block, respectively, and (b) the shear stress acting 
on the bottom surface of the rubber block, as a function 
of the distance xt the top surface has moved. Note that 
the bottom surface of the block is effectively pinned until 
the shear stress reaches a s « 0.44 MPa, corresponding 
to the friction coefficient /i = <t s /<to w 1.1. After the de- 
pinning, the velocity of the bottom surface of the block 
increases towards the velocity of the top surface, while 
the shear stress approaches a constant value determined 
by the kinetic friction coefficient. The physical reason for 
the peak in the shear stress at Xt ~ 0.07 cm is due to the 
flash temperature; the full flash temperature is not built 
up until the slip distance is of order the diameter of the 
macro asperity contact regions, i.e., of order 0.4 cm in 
the present case. 

Next, let us consider a case when the upper surface of 
the rubber block first accelerates with a = 500 m/s 2 for 
< t < 0.004 s, and then retards with a = -500 m/s 2 
for 0.004 < t < 0.008 s. In Fig. EJa) we show the 
shear stress, divided by the nominal pressure <7q, acting 
on the bottom surface of the rubber block, as a function 
of velocity vt of the top surface of the block. In Fig. 
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FIG. 19: A rubber block with thickness 0.5 cm sliding on an 
asphalt road. The nominal pressure ao = 0.4 MPa. The upper 
surface of the rubber block is clamped and undergoes uniform 
acceleration a = 500 m/s 2 for 0.004 s. (a) The velocities Vt 
and Vb of the top and bottom surface of the rubber block, 
respectively, as a function of the distance Xt the top surface 
has moved, (b) The shear stress acting on the bottom surface 
of the rubber block, as a function of x t . In the calculation 
the block is "divided" into N = 12 layers. 

I20f b) we show the velocity i>b of the bottom surface of 
the rubber block, as a function of time. 

9 Comparison with experiment 

A lot of experimental data has been presented in the 
literature related to rubber friction. However, in or- 
der to quantitatively compare the rubber friction the- 
ory with experimental data, both the viscoelastic mod- 
ulus E(ui,T) and the substrate surface roughness power 
spectrum C(q) must be known, but this information was 
never reported in any experimental investigations of rub- 
ber friction which I am aware of. Thus, in the present 
section we will only consider general (universal) aspects 
of rubber friction, and show that the theory is in good 
qualitative agreement with the presented experimental 
data. 

Let us first consider stationary sliding. Fig. 1211 shows 
the measured kinetic rubber friction coefficient as a func- 
tion of the logarithm of the sliding velocity, for a tread 
rubber block sliding on an asphalt surface 0. The veloc- 
ity dependence of Hk(v) is in good qualitative agreement 
with the theory (compare Fig. [2]with Fig. Illf a')). More 
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FIG. 20: A rubber block with thickness 0.5 cm sliding on an 
asphalt road. The nominal pressure ao = 0.4 MPa. The upper 
surface of the rubber block is clamped and undergoes uniform 
acceleration a = 500 m/s 2 for t < 0.004 s, and then uniform 
retardation a = -500 m/s 2 for 0.004 < t < 0.008 s. (a) The 
shear stress, divided by the nominal pressure ao, acting on the 
bottom surface of the rubber block, as a function of velocity 
v t of the top surface of the block, (b) The velocity Vb of the 
bottom surface of the rubber block, as a function of time. In 
the calculation the block is "divided" into N — 12 layers. 



generally, I have found that the theory, and all the exper- 
iments known to me, gives a maximal friction coefficient 
of order unity, and the position of the maximum in the 
range 10 -3 — 1CP 1 m/s, depending on the rubber com- 
pound and the substrate surface. In the absence of the 
flash temperature, according to the theory the maximum 
would instead occur at much higher sliding velocities, 
typically in the range 10 2 - 10 4 m/s. This illustrates the 
crucial role of the flash temperature. 

Next, let us consider the dependence of the friction 
coefficient fi on the background temperature. Fig. l2*2l 
shows exp erimental data (obtained using a portable skid 
tester [l^, [lj]) for a tread rubber block sliding on a road 
surface. The theory data in the same figure is calcu- 
lated for a rubber block sliding on an asphalt road sur- 
face (with the power spectrum given in Fig. [2J surface 
1) at v = 1 m/s. The viscoelastic modulus of the rubber 
used in the calculation is for a tread rubber. As expected, 
the rubber friction decreases with increasing background 
temperature. 
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FIG. 21: The friction coefficient /jasa function of the loga- 
rithm of the sliding velocity. For a rubber tread block sliding 
on an asphalt road surface at T = 18 °C. (Courtesy of Olaf 
Lahayne). 
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FIG. 22: The friction coefficient fi as a function of the back- 
ground temperature, during sliding of a rubber block on a 
rough substrate. The experimental data was obtained using 
a portable skid tester (la.ll7| for tread rubber on road surface. 
The theory data was calculated using the viscoelastic mod- 
ulus of a standard tread rubber, sliding on an asphalt road 
surface (with the power spectrum given in Fig. surface 1) 
at v — 1 m/s. 
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FIG. 23: A soft rubber block (thickness 1.1 cm) sliding on a 
wet concrete surface at the nominal pressure 0.2 MPa. The 
variation of the effective friction coefficient with the velocity 
v of the drive is shown for non-uniform sliding (see inset) 
involving acceleration (a ps 5 m/s 2 ) followed by retardation 
(a w —5 m/s 2 ). From Ref. 19] (with permission). 



calculations 1241 has nothing to do with the static friction 
coefficient, but rather is a kinetic effect related to the 
finite sliding distance necessary in order to fully build 
up the flash temperature. In fact, rubber on rough sub- 
strates does not exhibit any static friction coefficient, but 
only a kinetic friction which (for small sliding velocities) 
decreases continuously with decreasing sliding velocity. 

Finally, let us check if the theory agrees with the ob- 
servation that tires for racer cars (or racer motorcycles) 
exhibit much larger friction (but also much larger wear) 
than tires for passenger cars. In Fig. I show the 
calculated kinetic friction coefficient, using the measured 
viscoelastic modulus of a rubber compound from racer 
tires, and for a passenger car tire, assuming all other 
conditions identical. The maximum kinetic friction is 
about 50% higher for the racer compound which is also 
what is typically found experimentally 20] . This calcula- 
tion indicates that even for racer tires the main origin of 
the friction is due to the internal damping of the rubber 
rather than an adhesive tire-road interaction. 



Let us now consider non-stationary sliding. Fig. 1251 
shows results for a (soft) rubber block (thickness 1.1 cm) 
sliding on a wet concrete surface at the nominal pressure 
0.2 MPa. The variation of the effective friction coefficient 
with the velocity v of the drive is shown for non-uniform 
sliding (see inset) involving acceleration (a ss 5 m/s 2 ) 
followed by retardation (a f=s — 5 m/s 2 )[l9|. Note that 
the experimental data is of the general form predicted 
by the theory (see Fig. I24f> with a "start-up" peak due 
to the flash-temperature; the rubber block must slide a 
distance of order the linear size of the macro asperity 
contact regions before the full flash temperature has been 
developed, and this is the origin of the "start-up" peak. 

Note that the "start-up" peak in Fig. [23 and in the 



10 On the origin of the short-distance cut-off 

The rubber friction theory presented above assumes 
that the friction is due to the internal friction of the 
rubber. That is, the substrate surface asperities exert 
forces on the rubber surface and result in pulsating de- 
formations of the rubber. The deformations cannot occur 
completely adiabatically, but result in transfer of energy 
into random heat motion in the rubber. In calculating 
this asperity-induced contribution to the friction I in- 
clude only the road surface roughness with wavevectors 
q < qi . Here I will discuss the origin of the short distance 
cut off gi. In principle, there may be several different ori- 
gins of gi. For example, if the road is covered by small 
particles, e.g., dust or sand particles, with typical diam- 
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FIG. 24: A rubber block (thickness 0.5cm) sliding on an as- 
phalt road surface at the pressure 0.2 MPa. The (calculated) 
variation of the effective friction coefficient with the velocity 
v t is shown for two different cases where the drive first ac- 
celerate (a « 5 m/s 2 and 50 m/s 2 ) followed by retardation 
(a » —5 m/s 2 and —50 m/s 2 ). 
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FIG. 25: The kinetic (steady state) rubber friction coefficient 
for typical tire tread compounds for racer car (top) and per- 
sonal car (bottom), sliding on an asphalt road surface. The 
background temperature T = 60 °C. 



eter D, then one may expect q\ ~ 1/D. Similarly, on a 
wet road, the water trapped in the surface cavities may 
act as an effective short-distance cut off|2l|. However, 
for clean dry road surfaces, I believe that the cut off may 
be determined by the rubber compound properties. 

If no short-distance cut off q\ would exist the flash 
temperature and the surface stresses in the contact areas 
would increase as we study the contact regions at higher 
and higher magnification. However, when the flash tem- 
perature and surface stresses becomes high enough, the 
rubber will degrade j2sj. In fact, measurements on bulk 
samples have shown that natural rubber rapidly ther- 
mally degrades already at T « 200 C. The surface layer 
of the tire rubber is exposed to oxygen and ozone and will 
degrade even faster than in the bulk, even if the tempera- 



ture would be the same. The thermal degradation results 
in a thin layer of modified rubber at the surface, and we 
will make the basic assumption that the deformation of 
the rubber on length scales shorter than the thickness 
of the modified layer, gives a negligible contribution to 
the tire-road friction. Thus, with respect to friction, wc 
assume that the modified surface layer acts as a "dead" 
layer. 

The formation of a modified surface layer is a ther- 
mally activated mechanic-chemical process involving the 
breaking of chemical bonds (and the formation of new 
bonds). The rate of bond-breaking (at the temperature 
T and the tensile stress a) is assumed to be given by the 
standard expression of activated processes: 

w « Woe -^(i—/-o)/k B T (n) 

where the attempt frequency wq is usually of order 
10 12 s _1 , and where AE is the activation energy for bond 
breaking, and a c the stress necessary for bond-breaking 
at zero temperature. The weakest (chemical) bonds in 
rubber cross-linked with sulfur, are multi-sulfur cross- 
links for which AE is of order ~ 2 — 3 eV. Assuming a 
typical cross-link density, and that the applied stress dis- 
tributes itself entirely on the cross-links (which requires 
so high temperature that the rubber is liquid-like between 
the cross-links), and assuming that the same force acts 
on all the crosslinks, one can estimate the stress a c to be 
of order ~ 1 GPa. However, in inhomogeneous materials 
such as rubber, there will be a large distribution of local 
forces acting on the cross-links, and the cross-links where 
the local stress is highest will in general break first. 

The contact time of a tire tread block in the footprint 
area is typically of order 0.01 s, and since the dead layer 
has been worn out in about ~ 100 contacts, the rate w 
must be of order 1 s _1 . Taking w ~ 1 s^ 1 gives 

AE(1 - a/a c )/k B T = ln(w /w) « 28 (72) 

Note that the RHS of this expression is not very sensitive 
to the value of w. We can use (72) as a criteria to deter- 
mining the optimum cut-off q\ : We determine q\ so that 
Eq. (72) is satisfied, where T is now the flash tempera- 
ture in the contact area and a the stress in the contact 
area. 

We have found that using AE w 1.25 eV and a c w 
100 MPa result in a friction coefficient in good agree- 
ment with experiment. Note that AE is smaller than 
the typical energy to break sulfur cross link (which is of 
order 2 — 3 eV) or the energy to break a C-C bond along 
a carbon chain (which is of order 3 — 4 eV). However, in 
the present case the bond breaking is likely to involve re- 
action with foreign molecules, e.g., ozone or oxygen, and 
the effective activation energy for such processes may be 
much smaller than when the bond break in vacuum. The 
value for a c is higher than the rupture stress of macro- 
scopic rubber blocks, but a c refer to the rupture stress 



of very small rubber volume elements, which is higher 
than for macroscopic rubber blocks due to the absence 
of (large) crack- like defects |23j. It is clear that the pro- 
cesses which determines the cut-off q\ are closely related 
to tire tread wear. 
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11 Summary and conclusion 

When a rubber block slips on a hard rough substrate, 
the substrate asperities will exert time-dependent defor- 
mations of the rubber surface resulting in viscoelastic 
energy dissipation in the rubber, which gives a contribu- 
tion to the sliding friction. Most surfaces of solids have 
roughness on many different length scales, and when cal- 
culating the friction force it is necessary to include the 
viscoelastic deformations on all length scales. The en- 
ergy dissipation will result in local heating of the rubber. 
Since the viscoelastic properties of rubber-like materials 
are extremely strongly temperature dependent, it is nec- 
essary to include the local temperature increase in the 
analysis. In this paper I have developed a theory which 
describe the influence of the flash temperature on rubber 
friction. At very low sliding velocity the temperature in- 
crease is negligible because of heat diffusion, but already 
for velocities of order 10~ 2 m/s the local heating may be 
very important, and I have shown that in a typical case 
the temperature increase result in a decrease in rubber 
friction with increasing sliding velocity for v > 0.01 m/s. 
This may result in stick-slip instabilities, and is of cru- 
cial importance in many practical applications, e.g., for 
the tire-road friction, and in particular for ABS-breaking 
systems. 
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Appendix A: Average size of a macro-asperity 
contact region 

Consider the contact between two solids with nomi- 
nally flat surfaces. Fig. [2H shows the contact area at 
increasing magnification (a)-(d). At low magnification 
£ < 1 it appears as if complete contact occur between 
the solids. The macro-asperity contact regions typi- 
cally appear for a magnification somewhere in the range 
Cm = Qm/lo ~ 2 — 5 depending on the substrate surface 
and the rubber compound, where A/Aq w 0.25 — 0.3 is 
slightly below the site percolation threshold (for a hexag- 
onal lattice). The macro-asperity contact regions (c) 
breaks up into smaller contact regions (d) as the mag- 
nification is increased. 

Let us now estimate the (average) size of the macro- 
asperity contact regions. Assume that the contact 
patches have the surface areas Aj (i = 1,...,N). Let 





FIG. 26: The contact area at increasing magnification (a)- 
(d). The macro-asperity contact area (c) breaks up into 
smaller contact areas (d) as the magnification is increased. 



us define the probability distribution 
PA = ^J2S(A-A i ) 

i 

where N is the number of contacting asperities. The 
average area of a contact patch 
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Assume that the distribution Ph of summit heights is a 
Gaussian, 



Ph = 



1 



-exp 



-h 2 /2h* 
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(2tt) 1 /2/ 1 * 

where h* is the root-mean-square amplitude of the sum- 
mit height fluctuations. Assume that we can neglect 
the (elastic) interaction between the macro-contact ar- 
eas. Thus, using the Hertz contact theory the (normal- 
ized) area of real contact is I^fij 
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imoRo I dh (h — d)Ph 

Id 



(A3) 



where Aq is the nominal contact area, Rq the radius of 
curvature of the asperity, no the number of asperities per 
unit area, and d the separation between the surfaces. The 
number of contacting asperities per unit area 
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Substituting (A2) in (A3) and (A4), and introducing the 
new integration variable £ = (h — d)/h* gives 
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and 
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where x = d/h* . Thus, the average macro-asperity con- 
tact area 
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We can estimate the concentration of asperities, no, the 
(average) radius of curvature R of the asperities, and the 
rms summit height fluctuation h* as follows. We expect 
the asperities to form a hexagonal-like (but somewhat 
disordered) distribution with the lattice constant A m = 
27r/g m so that 



n » 2/(AiV3) = 



2ttV3 



0.029^ 



The height profile along some axis x in the surface plane 
will oscillate with the wavelength of order w A m roughly 
as h(x) w /iiCOs(q m x) where /ii = ^/2/iq (where ft-o is the 
rms roughness amplitude). Thus, 

- » /l"(0) = = 

Finally, we expect the rms of the fluctuation in the sum- 
mit height to be somewhat smaller than ho- In what 
follows we use|H n = 0.023<2^, l/R a = 0.92/i g^, and 
h* = 0.63 ho. Using these results and defining a = ttR 2 
we obtain from (A5) and (A7): 
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(A8) 
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In Fig. |23 we show the radius R of a macro-asperity 
contact region as a function of the relative contact area 
AA/Aq. The numerical data are well approximated by 



Rq m = a + b(AA/A y 



where a = 0.526, b = 3.636 and c = 0.729. This fit- 
function is also shown in Fig. [57| In a typical case when 
a tread block is slipping on a road surface, A(£ m )/Ao = 
AA/Aq w 0.25. Using Fig. |3J|this gives the diameter of 
the macro- asperity contact regions 2R w 3/go- If go = 
600 m _1 this gives 2i? w 0.5 cm. 



[1] D. F. Moore, The Friction and Lubrication of Elastomers 
(Oxford: Pergamon, 1972) 




FIG. 27: The (average) radius 7? of a macro-asperity contact 
area as a function of the (normalized) projected contact area 
AA/Aq. 



[2] 
[3] 
[4] 



ir 

[8 

[9 
[10 

[11 

[12 

[13 

[14 
[15 



(A10) [16 



[17 
[18 
[19 



[20 
[21 

[22 
[23 



B. N. J. Persson, Surf. Sci. 401, 445 (1998). 
B.N.J. Persson, Phys. Rev. Lett. 89, 245502 (2002). 

A. Galliano, S. Bistac, and J. Schultz, J. Colloid In- 
terf. Sci. 265, 372 (2003); A. Galliano, S. Bistac and 
J. Schultz, J. Adhesion 79, 973 (2003). 

B. N.J. Persson, Eur. Phys. J. E8, 385 (2002). 

K. A. Grosch, in The Physics of Tire Traction: Theory 
and Experiment, D. F. Hays, A. L. Browne, Eds. (Plenum 
press, New York-London 1974), page 143. 
B.N.J. Persson and A.I. Volokitin, Phys. Rev. B65, 
134106 (2002). 

M. Kliippel, G. Heinrich, Rubber Chem. Technol. 73, 578 
(2000). 

B. N.J. Persson, J. Chem. Phys. 115, 3840 (2001). 

C. Yang, U. Tartaglino and B.N.J. Persson, Eur. Phys. 
J. E19, 47 (2006). 

S. Westermann, F. Petry, R. Boes and G. Thielen Tire 

Technology International - Annual Review 2004, P- 68 

Edited by UKIP Media & Events, UK. 

S. Westermann, F. Petry, R. Boes and G. Thielen, 

Kautschuk Gummi Kunststoffe 57, 645 (2004). 

M.L. Williams, R.F. Landel and J.D. Ferry, J. Am. 

Chem. Soc. 77, 3701 (1955). 

A. R. Payne, J. Appl. Polymer Sci. 16, 1191 (1972). 

G. Heinrich and M. Kliippel, Advances in Polymer Sci- 
ence 160, 1 (2002). 

R. Ohhara, International Polymer Science and Technol- 
ogy 23, 25 (1996). 

Y. Sakai, Jidosha Gijutsu 21, 640 (1967). 

Olaf Lahayne, private communication. 

T. Huemer, W.N. Liu, J. Eberhardsteiner, H.A. Mang 

and G. Meschke, KGK Kautschuk Gummi Kunststoffe 

54, 458 (2001). 

Paul Haney, The Racing & High-Performance Tire, Co- 
published by TV MOTERSPORTS and SAE (2003). 

B. N.J. Persson, U. Tartaglino, O. Albohr and E. Tosatti, 
Nat. Mater. 3, 882 (2004). 

B.N.J. Persson, O. Albohr, G. Heinrich and H. Ueba, J. 
Phys. Condens. Matter 17, R1071 (2005). 
A rough estimate of a c may be obtained from the en- 
ergy to propagate a crack very slowly in rubber. In 



this limiting case no viscoelastic energy dissipation takes 
place in front of the crack tip, and the cra ck p ropaga- 
tion energy (per unit area)(see, e.g., Ref. |2^P Go ~ 
2na a 2 c /E{0) giving a c w [G E(0)/2na ] 1/2 , where £(0) 
is the zero frequency elastic modulus and ao the crack tip 
radius. In a typical case for styrene butadiene copolymer 
Go » 30 MPa [see, e.g., A.N. Gent, Langmuir 12, 4492 
(1996)] and with E(0) w 2 MPa and a » 1 nm we get 
cr c w 100 MPa. 



[24] J. A. Greenwood, in Fundamentals of Friction, Macro- 
scopic and Microscopic Processes, Ed. by I.L. Singer and 
H.M. Pollack (Kluver, Dordrecht, 1992). 

[25] J. A. Greenwood and J.B.P. Williamson, Proc. Roy. Soc. 
A295, 300 (1966). 

[26] P.R. Nayak, ASME J. Lubrication Technology 93, 398 
(1971). 



